Constraints on Enhanced Dark Matter Annihilation from IceCube Results 
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Excesses on positron and electron fluxes measured by ATIC, and the PAMELA and Fermi-LAT 
telescopes can be explained by dark matter annihilation in our Galaxy. However, this requires large 
boosts on the dark matter annihilation rate. There are many possible enhancement mechanisms, 
such as the Sommcrfeld effect or the existence of dark matter clumps in our halo. If enhancements 
on the dark matter annihilation cross section are taking place, the dark matter annihilation in the 
core of the Earth should also be enhanced. Here we use recent results from the IceCube 40-string 
configuration to probe generic enhancement scenarios. We present results as a function of the dark 
matter-proton interaction cross section, a xp weighted by the branching fraction into neutrinos, f u „, 
as a function of a generic boost factor, Bf, which parametrizes the expected enhancement of the 
annihilation rate. We find that dark matter models which require annihilation enhancements of 
0(100) or more and that annihilate significantly into neutrinos are excluded as the explanation for 
these excesses. We also determine the boost range that can be probed by the full IceCube telescope. 

PACS numbers: 95.35. +d, 95.30.Cq, 11.30.Pb 
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I. INTRODUCTION 

Several recent observations in our galaxy have found an 
excess on the electron and/or positron fluxes pHH] . How- 
ever the origin of these events is yet not clear. If they 
are produced accordingly to standard physics, pulsars [1] 
might be their sources as well as they could be secondary 
particles created by shock accelerated hadrons [5]. An- 
other exciting possibility is that they can be a signature 
of new physics. It was shown that dark matter annihi- 
lation in the galaxy can describe the data [6]. However 
most of these models require an enhancement on its anni- 
hilation rate. An enhancement mechanism can be found 
from a variety of possible phenomena, such as the Som- 
merfeld effect, the existence of dark matter substructures 
in the galactic halo, or a combination of both [7J [5] - 

Constraints on a boost factor on the dark matter anni- 
hilation rate have been recently derived from the Fermi- 
LAT diffuse gamma-ray measurement [S] and from their 
analysis of Milky Way satellite galaxies [28] . An indepen- 
dent and neat analysis has been performed by the authors 
of Ref [TO] , using the fact that dark matter annihilations 
at recombination time, redshift ~ 1000, would have in- 
jected secondary particles that would have affected the 
recombination processes. The measured power spectrum 
of the cosmic microwave background (CMB) can thus be 
used to set limits on the strength of the dark matter self- 
annihilation cross section. The IceCube collaboration has 
also performed an analysis searching for dark matter an- 
nihilations in the galactic halo by searching from an ex- 
cess diffuse neutrino flux over the expected atmospheric 
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flux [TT] . The results can also be used to set limits on a 
boost in the annihilation cross section. 

These methods show some degree of complementarity 
since each one provides better sensitivity to a different 
range of the mass of the dark matter particles or to dif- 
ferent annihilation channels. They also rest on different 
assumptions and approximations. Fermi and the CMB 
method are competitive in the low mass range (dark mat- 
ter massed < 1 TeV), while IceCube reaches 10 TeV, 
and can also probe annihilation directly into neutrinos. 
There is some halo-model dependency when setting lim- 
its on the annihilation cross section from the observations 
of satellite galaxies: the expected signal depends on the 
degree of cuspiness of the assumed halo. On the other 
hand, the CMB analysis does not depend on the shape 
of dark matter haloes, since there are no gravitationally 
bound structures at the redshift considered, but it de- 
pends on assumptions about the fraction of the energy 
released by the annihilating dark matter particles and 
how it is absorbed by the surrounding medium. 

In this paper we follow the alternative approach pro- 
posed in Ref. [J2]. The authors argue that an enhance- 
ment of the dark matter annihilation rate should also 
boost the neutrino flux from dark matter annihilations 
in the center of the Earth. A thermal annihilation cross 
section implies that dark matter capture and annihilation 
has not yet reached equilibrium in the Earth. However 
it is not necessary that the annihilation cross section has 
to remain the same as during the freeze-out period. If 
the annihilation rate is somehow enhanced in the post 
freeze-out period, the equilibrium might already have 
been achieved. As it will be shown in the next Section, 
the annihilation rate reaches its maximum at the equilib- 
rium, and then depends only in the capture rate. A boost 
in the annihilation rate increases the flux of annihilation 
products until it reaches its maximum. In this case the 
neutrino flux from the center of the Earth will be large 



enough to be detected by telescopes such as IceCube. In 
short, if the dark matter annihilation rate is enhanced, 
the timescale for equilibrium diminishes and the flux of 
annihilation products is expected to be much larger than 
the one away from equilibrium. 

We use the recently published IceCube results on a 
search for a diffuse flux of muon neutrinos [13] to set lim- 
its on a boost factor in the dark matter annihilation cross 
section. IceCube has measured neutrinos coming from 
near or below the horizon in the energy range 332 GeV 
and 84 TeV using data taken with the 40-string detec- 
tor configuration (IceCube-40). The analysis includes all 
events coming from near or below the horizon, and the 
result is compatible with the expected atmospheric neu- 
trino flux (see also Ref. [14 ). It is also generic enough to 
allow comparison with our prediction of the flux of muon 
neutrinos produced in dark matter annihilations in the 
center of the Earth. We determine this flux by simulating 
the annihilation of WIMP-type particles in the center of 
the Earth and propagating the neutrinos to the detector. 
A significant neutrino flux from dark matter annihilation 
should be seen above the expected atmospheric neutrino 
flux. If not, limits can be set on the model used to de- 
termine the dark matter signal. Our analysis shows that 
models which require very large boosts on the dark mat- 
ter annihilation rate in order to explain the excess seen 
in the galactic positron and electron flux, and have an 
annihilation channel into neutrinos, are ruled out. 

In the next Section we describe dark matter capture 
and annihilation in the Earth. In Sections Mil and IIVI we 
describe the signal simulation and calculate the expected 
number of events in the IceCube-40 detector. We then 
compare our results with the IceCube-40 published re- 
sults, showing the boost factor range that is excluded. 
Finally, in Section [V] we estimate the sensitivity region 
for the full 86-string detector. 



II. DARK MATTER ANNIHILATION IN THE 
EARTH 

Dark matter interactions in the Earth will be dom- 
inated by spin-independent elastic scattering, since the 
most abundant isotopes of the Earth core and mantle are 
spin nuclei. The time evolution of the number of dark 
matter particles will result from a balance between the 
capture (re) and annihilation rate {T A )- 

N = T c - 2T A . (1) 

The Earth dark matter capture rate is given by [T2] [15] : 
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here a xp is the spin-independent dark matter cross sec- 
tion off protons, v x and p x are the dark matter velocity 
and density in the halo and m x is the dark matter par- 
ticle mass. In deriving the above equation it has been 
assumed that for masses m x much larger than the nu- 
cleus mass, the dark matter interaction cross section off 
a proton is the same as off a neutron. In this case, the 
spin-independent cross section off a nucleus with mass 
number A is given by <j x n ~ A 4 a xp {l — 2m p /m x ), where 
m p is the proton mass. 

The annihilation rate depends both on the relative 
velocity scaled cross section (a A v) as well as on the 
dark matter distribution in the Earth [16] . The lat- 
ter can be given in terms of the parameter C A = 
(a A v)/V e ff, where T A = N 2 C A /2, and V e}f = 5.7 x 
10 22 (TeV/m x ) 3 / 2 cm 3 represents the dark matter effec- 
tive volume in the core of the Earth, assuming an isother- 
mal distribution [HI [17] . 

The solution to the dark matter time evolution in the 
Earth (Equation [I]) is then, 



T a = — tanh I — 
2 T 



(3) 



where £® is the age of the solar system and r = l/\/T c G A 
is the timescale for equilibrium between capture and an- 
nihilation. An enhancement on the annihilation rate will 
only be effective if it can accelerate equilibrium within the 
Earth. When the equilibrium stage is reached, the anni- 
hilation is maximum and depends entirely on the capture 
rate (r^ = Tc/2). Thermal relic dark matter candi- 
dates typically have (<J A v) r = 3 x 10~ 26 cm 3 s _1 , which 
makes C Ar ~ 5.3 x 10- 49 (m x /TeV) 3 ^ 2 s^ 1 and for these 
values the Earth is today far from equilibrium. 

In our analysis we consider scenarios where the annihi- 
lation cross section is enhanced by boosting the thermal 
relic annihilation rate T Ar by a generic factor Bf which 
affects the non-equilibrium rate through C A — BfC Ar . 
Such parametrization, though adequate to probe en- 
hancements due to Sommerfcld effect or to new inter- 
action mechanisms, cannot probe enhancements due to 
a possible dark matter halo substructure. In this latter 
case, a standard annihilation cross section could well ac- 
count for any possible signal, which would be just due to 
the increased local dark matter density, and not due to 
any new feature of the annihilation process itself. 



III. 



NEUTRINO PRODUCTION 
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The escape velocity at the Earth core is v e « 15 km/s 
and dark matter moves very slowly inside the Earth. 
The annihilation products will therefore be monochro- 
matic, and produced with the same energy as the dark 
matter mass. Here we consider the annihilation into a 
muon neutrino pair (which we call primary neutrinos). 
"Secondary" neutrinos are also produced in the decay 
of other primary annihilation particles, such as r's, t's, 



W's and b's. Since the explanation of the observed lep- 
tonic excess in terms of dark matter has also to ac- 
count for the fact that no antiproton excess was found 
by PAMELA [T5], annihilation into leptons is preferred 
when compared to hadrons. Secondary neutrinos from 
annihilation into charged lepton states, specifically on 
tt were analyzed in Ref. [12] . The energies of these sec- 
ondary neutrinos will be spread at relatively low values 
(~ 50 GeV) compared to the primary neutrino flux, and 
detection in neutrino telescopes is then disfavored (un- 
reasonably large boost factors would be needed to bring 
such flux over the atmospheric neutrino background to 
a detectable level). We will therefore not take into ac- 
count the secondary neutrino flux in our calculations, and 
concentrate in the easily detectable monochromatic flux 
from direct annihilations. 

The primary neutrino flux, produced from dark matter 
annihilations in the Earth's center is given by: 



M v = 500 GeV 
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dE v dAdtdfl ~ 4 tt Rl dE v 
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(4) 



where f v p is the annihilation branching ratio into v^V^, 
i?© is the radius of the Earth, and dN v /dE v is the en- 
ergy distribution of the neutrinos produced in the an- 
nihilations. We show our results for two generic cases, 
m x = 500 GeV and m x = 1000 GeV, which are repre- 
sentative of the models which fit the observed positron 
and electron excess, and later on for m x = 5000 GeV 
when comparing our results to others. Since at these en- 
ergies neutrinos practically do not lose energy in their 
way from the center of the Earth to the detector, the 
term dN v /dE v will then be a delta function at the dark 
matter masses considered. 



IV. CONSTRAINTS ON THE ANNIHILATION 
BOOST FACTOR FROM ICECUBE-40 



In order to predict the muon flux from dark matter 
annihilation in the Earth at the IceCube detector in 
the South Pole, we use the publicly available WimpSim 
code |19) . We simulate a monochromatic muon neutrino 
beam at the center of the Earth, with energy equal to the 
dark matter mass, by selecting the muon neutrino chan- 
nel. WimpSim simulates the propagation of these neu- 
trinos, including energy losses and charged and neutral 
current interactions as well as oscillations through the 
Earth. The muon neutrino flux at the detector is given 
as an output. Although we have chosen to simulate the 
flux specifically for the location of the IceCube detector, 
the detector site is not relevant in this case, and the re- 
sults and sensitivities presented in the next sections can 
easily be interpreted for any neutrino telescope of similar 
size. 

The number of muons from dark matter annihilation 
from a given angular region Q and during an exposure 




FIG. 1. Angular distributions of the primary muon neu- 
trino and antineutrino flux from dark matter annihilations at 
the center of the Earth with respect to the vertical direction 
(pink lines). The green curves show the angular distribution 
smeared by the detector angular resolution, taken as 2°. The 
bottom plot shows the distribution for 1 TeV dark matter, 
while the plot on the top for 500 GeV dark matter. Full lines 
are for neutrinos and dashed are for antineutrinos. 



t cxp in IceCube-40, is then: 



N u 



dE„dAdtdfl 



dE u t e xp A c jf 1 1. 



(5) 



We proceed then by convoluting the muon neutrino 
flux with the IceCube-40 effective area published in 
Ref. [13] . The effective area accounts for the detector ef- 
ficiency including the neutrino-nucleon interaction prob- 
ability, the muon energy loss from its production point 
to the detector, and the detector trigger and analysis ef- 
ficiency. As the neutrino beam is monochromatic, we 
use the corresponding value of A e ff for each dark matter 
mass considered. The muon-neutrino angular distribu- 
tion is the main parameter for background reduction in 
our analysis. Figure [T] shows the angular distribution 
from neutrinos from dark matter annihilation in the cen- 
ter of the Earth. It is collimated in a angle of less than 
approximately 3° around the vertical direction, 9 = 180° . 
This figure also shows the smearing effect due to the an- 
gular resolution of the detector, taken here conservatively 
to be 2° for up-going vertical events. We use therefore 
the effective area for the zenith range 150° < 6 < 180° 
from Ref. [13] . Even though such angular range is much 
wider than that from the expected signal, the effective 
area does not vary considerably with zenith angle below 
about 10 TeV, where Earth absorption effects start to be 
important. 

For our analysis, we choose an angular window from 
the vertical direction that contains 90% of the expected 
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FIG. 2. Predicted number of v^ + v^ from dark matter an- 
nihilation in the Earth at IceCube-40. Results are for 375.5 
days of exposure for different boost factor values (color coded 
as labeled). Dashed lines are for 1 TeV and solid lines for 500 
GeV dark matter. 



signal (4.1°(3.7°) for 500(1000) GeV dark matter respec- 
tively, where the angular experimental resolution has 
been taken into account), and use the calculated num- 
ber of events with Eq. [5] as signal. We have not used any 
energy information in this analysis, just using the total 
number of events predicted and detected to derive our ex- 
clusion regions. We will however present a study on the 
effect of the detector energy resolution in the sensitivity 
study for the complete IceCube detector in the next sec- 
tion. Figure [2] shows the predicted number of v^ + V p 
in IceCube-40, with t cxp = 375.5 days as a function of 
the dark matter nucleon cross section a xp scaled by the 
branching ratio for different boosts factors. As can be 
seen, larger enhancements on the annihilation rate re- 
sult in a shorter time for equilibrium among capture and 
annihilation. 

In order to study the model rejection power of our anal- 
ysis in a quantitative way, we determine the statistical 
significance of the predicted dark matter signal, S/VB~, 
as a function of boost factor, where S is the number of 
signal events predicted in the given angular cone. As 
background, B, we use the measured number of events 
in IceCube-40 in the same angular regions. 

Two recent IceCube publications, a measurement of 
the atmospheric neutrino flux |14j and a search for a dif- 
fuse E~ 2 flux of cosmic origin [13], give results that are 
consistent with the expected atmospheric neutrino flux. 
We use the publicly available data from the diffuse anal- 
ysis which is available at Ref. [20]. This analysis selects 
events coming from near or below the horizon and their 
background is composed by atmospheric muons arriving 
from above the detector and misreconstructed as an event 



coming from below. The rejection of these events is de- 
scribed in Ref. |13j and the background contamination 
in the data sample is estimated to be less than 1%. The 
final data sample contains 13K events and it is given as 
a function of the zenith angle. This allows us to select 
events which come in the same direction as expected sig- 
nal dark matter annihilations. The data is then reduced 
to 14 (9) events when considering the angular regions ex- 
pected for 500 (1000) GeV dark matter neutrinos. We 
compare this number to the predicted number of events 
from our signal choices and for different boost factors. 

The result is presented in Fig.[3j which shows the boost 
factor versus a xp scaled by the vv branching ratio region 
which was probed in this analysis. The plot to the left 
is for 500 GeV dark matter and the one to the right is 
for 1 TeV. The light green full line shows a 5a exclusion 
and above this line the exclusion is more strict according 
to the color coded bar on the right side. This bar indi- 
cates the level of statistical significance. The best limit 



on <7 VT) comes from the XenonlOO collaboration 



At 



90% CL 500 GeV dark matter models are constrained 
above a xv « 4 x 10 -44 cm 2 and 1 TeV dark matter above 
a xp « 8 x 10 -44 cm 2 . We draw this limit on Fig.wHust as 
a reference, since direct detection results do not depend 
on / y j7 neither on any boost factor. For such limit, boost 
factors above 215 and 58 respectively are excluded at a 
5cr level. This exclusion requires, however, that fuv = 1, 
ie., annihilation exclusively into neutrinos, which may 
be difficult to justify phenomenologically. But it should 
be noted that the excluded region in Fig. [3] does not re- 
quire necessarily a 100% annihilation into neutrinos. One 
can lower the branching factor by increasing the required 
boost factor limit. As one moves down along the green 5<r 
limit line, larger boost factors are disfavored, as the quan- 
tity a xp x f vV becomes lower. Lower values of a xv x f vV 
can be achieved either by models with relatively high 
cross sections but a lower branching ratio to neutrinos, 
or by lower cross sections but a higher annihilation prob- 
ability to neutrinos. In either case, the expected signal in 
IceCube would decrease and a higher boost factor would 
be needed to bring it to the current sensitivity of the de- 
tector. Higher boost factors to the right of the green line 
would then be disfavored at a higher significance. 

The degeneracy of the product o xp x j vV can be bro- 
ken when considering specific dark matter models with 
known branching ratio to neutrinos and cross section 
with protons. Fig. [3] can then be used to determine the 
minimum boost factor which is disfavored at a 5er level. 



V. PREDICTION FOR FULL ICECUBE 
DETECTOR 

The full IceCube detector with 86 strings (IceCube- 
86) is now completed and taking data since May 2011. 
In order to estimate the full detector sensitivity, we can 
use the fact that the atmospheric neutrino flux measured 
by IceCube-40 [131 [H] is consistent with model expecta- 



M^ = 500 GeV 



1000 GeV 




FIG. 3. Boost factor versus dark matter nucleon interaction cross section a xp scaled by vv branching ratio {fun)- The light 
green full line indicates the exclusion at ha level based on our analysis of IceCube-40 results. Above this line the exclusion 
is more strict accordingly to the color coded bar at the side. The color coded bar represents the statistical significance level 
(S/\/~B). Both the bar as well as the contour in this figure starts at 3<r. The black (red, blue) dashed line is the 5<r sensitivity 
for the IceCube-86 telescope with 50 (30, 10)% energy resolution. As a reference, the black full lines represent the 90% CL 
constraints on a xp from direct detection [23] (independent of the boost factor and of the annihilation branching ratio f v y). The 
plot to the left is for 500 GeV dark matter and the one to the right for 1000 GeV. 



tions as, for example, the one proposed by Ref. [3T]. We 
will also go one step further than in the calculations of 
the previous Section, and assume that IceCube will be 
able to estimate the neutrino energy with a given resolu- 
tion cte- We assume three benchmark energy resolutions, 
<je = 0.1E,0.3E and 0.5E which go from the very op- 
timistic to the more conservative situation in a neutrino 
telescope. The 40-string configuration has its energy res- 
olution between 50% and 80%. We assume here that the 
complete detector can reach a better resolution in any 
case. We note that energy estimation in neutrino tele- 
scopes is a difficult task since muon tracks above a few 
hunderd GeV will cross the detector volume, and the 
neutrino energy can only be estimated through model- 
dependent deconvolution methods. 

An effective area at trigger level for the complete 
IceCube-86 detector has been published in Ref. [25"] . 
However the effective area at final analysis level can dif- 
fer significantly from trigger level, since data quality cuts 
are applied to the data sample to reduce background, but 
also inevitably reducing signal efficiency. In practice, one 
can just rescale the IceCube-40 effective area by 2.15 
(which factors 40 to 86 strings), since at the energies 
considered in our analysis, and for vertical events, each 
string is practically an independent detector and there- 
fore, in a first approximation, the effective area scales 
proportionally to just the number of strings. We expect, 
though, dedicated analysis with IceCube-86 to be more 
efficient at lower energies than the IceCube-40 analysis 
we have used in the previous Section and, since we will 
be integrating A e ff x flux in a range of energies, we need 
to be careful with the behavior of the effective area with 
energy. We have therefore normalized the IceCube-86 ef- 
fective area from Ref. [S3] with the rescaled IceCube-40 



area at 10 TeV, where the IceCube-40 analysis is opti- 
mal. The shape of the IceCube-86 effective area at lower 
energies takes automatically into account the improved 
capabilities of the full detector at energies below 1 TeV. 
We have then all the ingredients we need to estimate 
the sensitivity of the IceCube-86 detector. We calculate 
the number of atmospheric neutrino events from within 
a vertical cone of 4.1° and 3.7° aperture by convoluting 
the detector effective area with the parametrization of 
the atmospheric neutrino flux taken from Ref. [21J . In- 
stead of choosing a delta function for the spectrum as in 
the previous case, we introduce a smearing in the energy 
according to the assumed energy resolution. In practice, 
this translates into that we perform the A e g x flux inte- 
gral between m x — ue and m x + ue- We assume that the 
measurement of IceCube will be compatible with the cal- 
culated number of background events and calculate the 
sensitivity to an excess neutrino flux from the center of 
the Earth under this assumption. The result of this exer- 
cise is shown in Fig. [3] The dashed black (red, blue) lines 
in both plots represent 5er sensitivity that can be achieved 
in one year live time of IceCube-86 (i C xp86 — 365 d) as- 
suming a 50 (30,10)% energy resolution. 



VI. 



CONCLUSIONS 



If dark matter annihilation in our galaxy is respon- 
sible for the excesses seen by ATIC 1], PAMELA [2 
and Fermi-LAT [3], an enhancement on the dark mat- 
ter annihilation rate in the Earth should also be fore- 
seen. This enhancement might bring the equilibrium 
among the dark matter capture and annihilation rates 
to occur much earlier than the timescale expected from 



a purely thermal annihilation rate. In this case the neu- 
trino flux from dark matter captured in the Earth can be 
quite large. We present results as a function of the dark 
matter-proton interaction cross section, a xp weighted by 
the branching fraction into neutrinos, f v y, as a function 
of a generic boost factor, Bp, which parametrizes the 
expected enhancement of the annihilation rate. In this 
sense, it is important to note that our results do not de- 
pend on the details of the mechanism which enhances the 
annihilation cross section. We have used two benchmark 
models, a 500 GeV and 1 TeV generic WIMP annihilat- 
ing in the center of the Earth to scan the (a xp x f vV , Bp) 
parameter space and set constrains on this 2-dimcnsional 
space using current IceCube results. Our calculations as- 
sume that the dark matter velocity distribution is Gaus- 
sian as well as that dark matter collected in the Earth's 
core follows an isothermal distribution. 

In order to explain the positron and electron excesses 
that are seen, models also have to cope with the fact that 
the antiproton spectrum measured by PAMELA [TB] is in 
full agreement with the expectation from secondary pro- 
duction of antiprotons from propagation of cosmic rays 
in the galaxy. This rules out as an explanation of the 
excess many dark matter models with preferred annihi- 
lation into heavy products, producing many antiprotons. 
Therefore leptophilic models, which propose dark mat- 
ter annihilation exclusively into leptons, are favored as 
an explanation to the excesses found. We have shown 
that, when 500 (1000) GeV dark matter annihilates into 
a large fraction of neutrinos, annihilation boost factors 
of the 0(100) and above are already constrained by our 
analysis at a 5<r level, or higher, depending on the interac- 
tion cross section assumed. Thus, leptophilic models [26] 
which favor primary neutrino production are constrained 
by our results. 

In order to compare this analysis to others, we show 
in Figure Rl limits on (a a X v) as a function of m x from 
Fermi [2STCMB [ID] and IceCube [25]. Our bounds as a 
function of (a a X v), can be determined from equations] 
It can also easily be visualized in Figure [3] where a 3<r 
significance is shown as the edge of the shaded area: the 
limit on the annihilation cross section is determined from 
each boost factor value on this curve. In Figure Rl these 
are shown for two choices of a xp , 3 x 10~ 44 crrr (blue 
stars) and 1 x 10 -44 cm 2 (blue solid squares), which are 
below the current Xenon limit [23]. 

Our results and IceCube's, are the only ones to 
probe annihilation into neutrinos. Other limits com- 
ing from searches with gammas from satellite galaxies by 
Fermi 28J or from the analysis based on the CMB [TDJ 
imprint at high redshifts, are complementary, since not 
only they probe different annihilation channels, but also 



rest on different underlying assumptions. 

We also investigated the reach of the completed Ice- 
Cube 86-string detector, and we present results of its 
sensitivity in the {a xp x f vV , Bp) parameter space. We 
have estimated how using neutrino energy information 
could improve the analysis, and in Fig. [3] we show the ex- 
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FIG. 4. Bounds on the dark matter annihilation cross section 
a a versus mass m x by Fermi [28], IceCube [29], analysis over 
CMB data [TU] and this work, where two different a xp values 
of 3 x 10" 44 cm 2 (blue star) and 1 x 10" 44 cm 2 (blue solid 
square) are assumed. Here we also plot our results assuming 
a 5 TeV dark matter mass. Our limits correspond to a 3a 
significance level, while the others are at 90% CL. 



pected sensitivity for three different assumptions of the 
detector energy resolution. 
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